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,—1 . ABSTRACT 

I ■ 

IL" ■ We report on the detection of very-high energy (VHE, E > 100 GeV) y-ray emission from NGC 1275, the central radio galaxy of the Perseus 
. cluster of galaxies. The source has been detected by the MAGIC telescopes with a statistical significance of 6.6 cr above 100 GeV in 46 hr 
k>( . of stereo observations carried out between August 2010 and February 2011. The measured differential energy spectrum between 70 GeV and 

■ 500 GeV can be described by a power law with a steep spectral index of F = -4.1 ± 0.7jMf + O.Sjvn, and the average flux above 100 GeV is 

■ = (1.3 ± Q.lstat ± 0.3jy„) X 10"" cm"^ s"' . These results, combined with the power-law spectrum measured in the first two years of observations 
■ " " ' by the Fermi-LAT above 100 MeV, with a spectral index of F -2.1, strongly suggest the presence of a break or cut-off around tens of GeV 

in the NGC 1275 spectrum. The light curve of the source above 100 GeV does not show hints of variability on a month time scale. Finally, we 
report on the nondetection in the present data of the radio galaxy IC 310, previously discovered by the Ferm;-LAT and MAGIC. The derived flux 



upper limit FJ; (> 300 GeV) = 1.2 x 10 cm s is a factor ~ 3 lower than the mean flux measured by MAGIC between October 2009 and 



7U.L., 

y 

February 2010, thus confirming the year time-scale variability of the source at VHE. 
Key words, galaxies: active — galaxies: jets — galaxies: individual (NGC 1275) — galaxies: individual (IC 310) — gamma rays: galaxies 



1. Introduction 

NGC 1275 (z = 0.0179), the central dominant galaxy of 
the Perseus cluster, harbors one of the closest active galac- 
tic nuclei ('AGN), already included in the original Seyfert 
list (ISevfertl Il943h . The AGN is a very bright radio source 
sho wing an extended jet with FanarofF-Rile y 1 m orphology 
(e.g. IVermeulen etall Il994t iButtiglione eTall l2010h . The op- 
tical emission of the nucleus is variable and strongly polarized 



from 3% to 6% ('Maza'.'l979l: lMartin et al.Lll983 h. implying that 
the relativistic jet contribute s sign ificantly to the optical con- 
tinuum f Angel & Stockman. 1980). The source has also been 
classified as a BL Lac object ( Veron , 1978,) . However, the jet in- 
creases its inclination from 10° to 20° on milliarcsecond scales 
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up to 40° to 60° at arcsecond scales (iKrichbaum et all Il992h . 
Due to its brightness and proximity this source is ideally suited 
to study the physics of relativistic outflows and t he "feedback" 
effbcts of the jet on th e cluster environment (e.g. iFabian et all 
l200ilGallagheai2009l) . 

In fact, NGC 1275 is one of the closest y-ray emitting AGN. 
It was first unambiguously detected in the high-energy (HE, 
100 MeV < E < 100 GeV) y-ray range by the Fermi Large 
Area Telescope (LAT) dAbdo et all l2009b . during the first four 
months of all-sky-survey observations, with an average flux 
above 100 MeV of Fy = (2.10 + 0.23) x lO"'' cm-^ s^^. The 
differential energy spectrum between 100 MeV and 25 GeV 
was described well by a power law with a spectral index of 
r - -2. 17 +0.05. While no variability was observed during these 
four months of observations, subsequent results based on the first 
year of Fermi-LAH observations ( Kataoka et al. , 20 1 0) show ev- 
idence of flux variability on time scales of months. Furthermore, 
the average y-ray spectrum show a significant deviation from a 
simple power law, indicating an exponential cut-off at the break 
photon energy of E, = (42.2 ±19.6) GeV. 

More recently, the resul ts obtained from the firs t two years 
of Fermi-LA^ observations (iBrown & Adamsll201 ih have given 
clear evidence for variability on time scales of days above 
800 MeV, revealing that several major flaring events occurred 
during the two-year observation period. A harder-when-brighter 
correlation between flux and spectral index was also found. 
Brighter and therefore harder > GeV states are then promis- 
ing for triggering observations at very high energy (VHE, E > 
100 GeV). Finally, pres e nt upp er limits at VH E provided by 
MAG IC-I (lAleksic et al.L 1201 Oah and VERITAS jAcciari et al.L 
[2()09| ) combined with the Ferm/-LAT results mentioned above 
suggested that NGC 1275 may have a break or cut-off in the 
spectrum around tens of GeV. 

The Perseus galaxy cluster contains another y-ray source, the 
radio galaxy IC 310, located at ~ 0.6° from NGC 1275. It was 
discovered in 2010 by the Fermi-LAT in the 30 GeV - 300 GeV 
energy range (Neronov et al., 2010) and for energies > 260 GeV 
by MAGIC (Aleksic et al. u2010b ). The combined MAGIC and 
Fermi-LAI spectrum is consistent with a flat spectral energy dis- 
tribution (SED) stretching without a break over more than three 
orders of magnitude in energy (2 GeV - 7 TeV). The spectrum at 
VHE measured by MAGIC has a spectral index of F = -2.00 + 
0.14, and the mean flux above 300 GeV, from October 2009 to 
February 2010, was Fy = (3.1 + 0.5) x 10"'- cm^^ g-i strong 
hints of a week to a year time-scale variabiUty were seen in the 
MAGIC data. 

In this letter we present the results of the observations of 
NGC 1275 at VHE performed with the MAGIC telescopes be- 
tween August 2010 and February 2011, which resulted in the 
first detection of the source above 100 GeV. The same ob- 
servational campaign also provided results on the variability 
at VHE of IC 310. This letter is accompanied by a separate 
paper dedicated to the study of the Perseus cluster environ- 
ment, foc using on possible VHE y-ray emission induced by cos- 
mic rays (lAleksic et al.L 1201 Ibl) . The multiwavelength emission 
of NGC 1275 from radio to VHE will be addressed in future 
work. 

2. Observations and analysis 

The MAGIC system consists of two 17 m dish Imaging Air 
Cherenkov Telescopes (lACTs) located at the Roque de los 
Muchachos observatory, in the Canary Island of La Palma 
(28.8°N, 17.8°W, 2200 m a.s.l.). Since late 2009 the telescopes 



are working in stereoscopic mode providing an excellent sen- 
sitivity of < 0.8% of the Crab flux (C.U.f] for energies above 
~ 300 GeV in 50 hr of observations (Aleksic et al., 201 la) and 
a trigger energy threshold of 50 GeV, which is the lowest among 
the existing lACTs. The MAGIC telescopes are currently the 
most sensitive instrument between 50 GeV and 200 GeV, allow- 
ing us to extend up to the TeV scale the observations carried out 
by the FermZ-LAT 

The Perseus galaxy cluster region was observed by the 
MAGIC telescopes during two campaigns. The first one was car- 
ried out between October 2009 and February 2010, for a total ob- 
servation time of 45.3 hr This survey result ed in the discovery o f 
the radio galaxy IC 310 as VHE emitter ( Aleksic et al.L l2010bl). 
The latest campaign (total observation time of 53.6 hr), which 
resulted in the detection of NGC 1275 at VHE presented in this 
letter, was performed between August 2010 and F ebruary 201 1 . 
The source was observed in the wobble mode (F omin et al.L 
1994), with data equally split in four pointing positions located 
symmetrically at 0.4° from NGC 1275, in order to ensure opti- 
mum sky coverage and background estimation. The survey was 
carried out during dark time at low zenith angles (from 12° to 
36°), which guaranteed the lowest energy threshold (~ 50 GeV). 

The data analysis was performed using t h e standard 
MAG IC software package (lAlbert et al.L l2008bl; lAUu et all 
'2009), taking advantage of newly developed stereoscopic 
analysis routine s (Moralejo et al., 2009; Aleksic et al., 201 lal 
iLombardi et al.L 1201 ll) . The analysis cuts applied to NGC 1275 
data were optimized by means of contemporaneous Crab Nebula 
data and Monte Carlo (MC) simulations. 

After the application of standard quality checks, 7.9 hr of 
data were rejected mainly due to nonoptimal atmospheric con- 
ditions. The selected sample used for deriving the results pre- 
sented here is therefore composed by 45.7 hr of good quality 
stereo data. 

3. Results 

The 0^ distribution^ with respect the signal region and the back- 
ground (estimated from 3 distinct regions), for energies above 
100 GeV, are shown in Fig. [1] We found an excess of 522 + 81 
events, corresponding to a significance of 6.6 standard devi ations 
(o-), calculated according to the Eq. 17 of iLi & Mai (Il983h . It is 
worth noting that the background estimation is not affected by 
a possible IC 310 y-ray contribution, since the latter source was 
not detected in the present data. 

The NGC 1275 differential energy spectrum measured by 
MAGIC between 70 GeV and 500 GeV can be described by a 
simple power law (x^/rtdof = 0.76/1) 

^=(3.1 + 1.0..±0.7.,„.„)xl0-(^^f, (1) 

in units of cm"^ s"' TeV"', with F = -4.1 + OJ^tat ± 
03sysE The mean flux above 100 GeV is Fy = (1.3 ± 

' In this letter C.U. stands for Crab u nits, defined a s the fra ction of 
the Crab Nebula flux given in Eq. 1 of lAleksic et"an d2011ah . which 
corresponds for energies above 100 GeV to 5.4 x 10"'" cm"^ s"'. 

- The 9^ is the squared angular distance betw een the arrival dire ction 
of the events and a given nominal position (e.g. lDaum et aliri997h . 

^ The systematic errors of the flux normalization and the energy 
spectral slope considered here have been estimated to be 23% and 
±0.3, respectively, whereas the systematic error on the energy scale 
is 17%. These values are more conservative than those presented 
m Aleksi c et al.l ( 1201 Ij) . given the flux weakness and the spectral steep- 
ness of NGC 1275, as measured by MAGIC. 
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Fig.l. distributions of the NGC 1275 signal and the back- 
ground estimation from 45.7 hr of MAGIC stereo observations 
taken between August 2010 and February 201 1, above an energy 
threshold of 100 GeV. The region between zero and the vertical 
dashed line (at 0.02 degrees^) represents the signal region. 



0.2j,a, + 0.3.5y.s;) X 10"" cm^^s"', corresponding to (2.5 + 
QAstai ± 0.6jyj,)% C.U. The steepness of the spectral index mea- 
sured by MAGIC strongly supports the presence of a break or 
cut-ofF in the NGC 1275 spectrum around tens of GeV, as al- 
ready suggested by the F ermi-LACT results (IKataoka et al.L 120101: 
iBrown & Adamsl 1201 ll) . and is consistent with the upper hmits 
on the flux at V HE provided by MAG IC-I (lAleksic et al.Ll2010al) 
and VERITAS (lAcciari et all l2009l) . The rapid decline in the 
spectrum, which causes the NGC 1275 signal to vanish above 
approximately 500 GeV, permits investigating possible VHE y- 
ray emissions induced by cos mic rays in the Perseu s cluster en- 
vironment above that energy (lAleksic etal.L 1201 Ibl) . 

In Fig. m the SED measured by MAGIC is com- 
pared with the results in the 100 MeV - 100 GeV 
range provided by the Fermi-LAT, ave raging Fermi data 
over t he first yea r (Kataoka et alt 1201 O l) and the first two 
years jBrown & A dams, 2011). The comparison suggests that 
a significant spectral steepening occurs around ~ 100 GeV. 
However, the present non-simultaneous data do not allow dis- 
cussing whether the spectral change corresponds to a break be- 
tween two power laws or exponential cut-off". 

The August 2010 to Februai-y 201 1 Hght curve of NGC 1275 
computed for an energy threshold of 100 GeV and with a 
monthly binning is shown in Fig. [3] No evidence of variability 
can be derived from these measurements. In fact, fitting the light 
curve with a constant flux hypothesis yields ax^/ndof - 7.4/6, 
corresponding to a probability P{x^) = 0.29. 

The significance skymap of the central region of the Perseus 
cluster above 100 GeV is shown in Fig. |4] A hot spot at 
> 6 cr significance level consistent with the sky position of 
NGC 1275 is present. The position of the radio galaxy IC 310 
is also shown. No significant excess events coming from IC 310 
have been found in the observations presented here. The corre- 
sponding integral flux upper limit above 300 GeV (performed 
using the iRolke et all 2005 method, with a confidence level 
of 95%, and a total systematic uncertainty of 30%) is Fy -^ {> 
300 GeV) = 1.2 x 10"'^ cm^^ s"', for a spectral index of 
r - -2.0 (i.e. the spectral index of the source previously mea- 
sured by MAGIC). This value is about a factor 3 lower than 
the average integral flux Fy{> 300 GeV) = (3.1 + 0.5) x 
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Fig. 2. NGC 1275 spectral energy distribution measured by 
MAGIC between 70 GeV and 500 GeV (filled squares), together 
with the resu lts above 100 MeV achieved from the first year 
(filled circles, Kataoka et al.','2010'), and from the first two years 
(open circles. Brown & Adams, 2011) of the Ferm/-LAT all- 
sky-survey observations. The power-law fits to the ferm/-LAT 
data (extrapolated up to 1 TeV) are also shown, together wit h 
the exponential power-law fit provided in lKataoka et alj (l2010l) . 
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Fig. 3. NGC 1275 Hght curve between August 2010 and 
February 201 1 above an energy threshold of 100 GeV, and with a 
month time-scale binning. No hints of variability are seen in the 
data. The dashed horizontal line represents the constant func- 
tion resulting from the fit to the data. For the December 2010, 
January 2011, and February 2011 data, the upper limits on the 
flux abov e 100 GeV for a spe ctral index of F = -4.0 (calculated 
using the iRoIke et all l2005l method with a confidence level of 
95%, and a total systematic uncertainty of 30%) are also shown 
(open dashed arrows). 

10 cm"^ s ' measured by MAGIC from October 2009 to 
February 2010 dAleksic et al.L 1201 Obi) , thereby confirming the 
variability of the latter source on a year's time scale. 

4. Conclusions 

The MAGIC telescopes have detected VHE y-ray emission from 
NGC 1275, the central radio galaxy in the Perseus cluster, at a 
statistical significance of 6.6 cr from observations performed be- 
tween August 2010 and February 2011. The corresponding av- 
erage flux above 100 GeV is Fy - (1.3 + 0.2^,0, + 0.3jv.s/) x 
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Fig. 4. Significance skymap of the central region of the Perseus 
galaxy cluster from 45.7 hr of MAGIC stereo observations taken 
between August 2010 and February 2011, above an energy 
threshold of 100 GeV. The NGC 1275 position is marked with 
a black cross, whereas the position of the radio galaxy IC 310 is 
shown by a green star The PSF of about 0.12° is also displayed. 
The hot spot in the map at a significance level > 6 cr is consis- 
tent with a point-like emission coming from the NGC 1275 sky 
position. 



10 cm ^ s '. This is the fourth nearby ra dio galaxy detected 
at VHE, after M 8 7 ( Aharo nian etdl I2003L 120 06: Albert et al., 
l2008aH Acciari et al.. 2008), Cen A (lAharonian et al.. ,20091 and 
IC 310^1eksic et al., 20109)- The MAGIC observation yields a 
spectrum that can be fitted between 70 GeV and 500 GeV by a 
simple power law with a spectral index of F = -4.1 + Q.l stat ± 
03svs t- This result, combined with previous Fermi-LKT re- 
sults (lAbdo et al.L l2009t iKataoka et al.1 l2010HBrown & Adamsl 
1201 Ih . showing a power-law spectrum with a spectral index of 
r ^ -2.1 above 100 MeV, strongly suggests the presence of a 
break or cut-off around tens of GeV in the NGC 1275 spectrum. 
No evidence of variability on month time scale has been found 
above 100 GeV. Finally, the variability on a year time scale of 
the source IC 310 (Neronov et al., 2010; Aleksic et al., 2010b) 
has been confirmed. The upper limit above 300 GeV presented 
here is in fact about a factor 3 lower than the flux measured by 
MAGIC between October 2009 and February 2010. 



Acknowledgements. We would like to thank the Instituto de Astrofi'sica de 
Canarias for the excellent working conditions at the Observatorio del Roque de 
los Muchachos in La Palma. The support of the German BMBF and MPG, the 
Italian INFN, the Swiss National Fund SNF, and the Spanish MICINN is grate- 
fully acknowledged. This work was also supported by the Marie Curie program, 
by the CPAN CSD2007-00042 and MultiDark CSD2009-00064 projects of the 
Spanish Consolider-Ingenio 2010 program, by grant DO02-353 of the Bulgarian 
NSF, by grant 127740 of the Academy of Finland, by the YIP of the Helmholtz 
Gemeinschaft, by the DFG Cluster of Excellence "Origin and Structure of 
the Universe", by the DFG Collaborative Research Centers SFB823/C4 and 
SFB876/C3, and by the Polish MNiSzW grant 745/N-HESS-MAGIC/2010/0. 
CP. gratefully acknowledges financial support of the Klaus Tschira Foundation. 



References 

Abdo, A. A., Ackermann, M., Ajello, M., et al. 2009, ApJ, 699, L31 
Acciari, V. A., Beilicke, M., Blaylock, G., et al. 2008, ApJ, 679, L397 
Acciari, V. A., Aliu, E., Arlen, T., et al. 2009, ApJ, 706, L275 
Aharonian, F. A., Akhperjanian, A. G., BeiUcke, M., et al. 2003, A&A, 403, LI 
Aharonian, F. A., Akhperjanian, A. G., Bazer-Bachi, A. R., et al. 2006, Science, 
314, 1424 

Aharonian, F. A., Akhperjanian, A. G., Anton, G., et al. 2009, ApJ, 695, L40 
Albert, J., AUu, E., Anderhub, H., et al. 2008a, ApJ, 685, L23 



Albert, J., Aliu, E., Anderhub, H., et al. 2008b, ApJ, 674, L1037 
Aleksic, J., Antonelli, L. A., Antoranz, P., et al. 2010a, ApJ, 710, L634 
Aleksic, J., Antonelli, L. A., Antoranz, P., et al. 2010b, ApJ, 723, L207 
Aleksic, J., Alvarez, E. A., Antonelli, L. A., et al. 2011a, ar Xiv:1108.1477t sub- 
mitted to Astroparticle Physics 
Aleksic, J., Alvarez, E. A., Antonelli, L. A., et al. 201 lb. iarXiv:1111.5544l sub- 
mitted to A&A 

AUu, E., Anderhub, H., AntonelU, L. A., et al. 2009, Astroparticle Physics, 30, 
293 

Angel, J. R. R, & Stockman, H. S. 1980, Ann. Rev. Astron. Astroph., 18, 321 
Brown, A. M., & Adams, J. 201 1, MNRAS, 413, 2785-2790 
Buttiglione, S., Capetti, A., Celotti, A., et al. 2010, A&A, 509, A6 
Daum, A., Hermann, G., HeB, M., et al. 1997, Astroparticle Physics, 8, 1 
Fabian, A. C, Johnstone, R. M., Sanders, J. S., et al. 2008, Nature, 454, 968-970 
Fomin, V. P., Stepanian, A. A., Lamb, R. C, et al. 1994, Astroparticle Physics, 
2, 137 

Gallagher, J. S., 2009, AN 330, 1040G 

Kataoka, J., Stawarz, L., Cheung, C. C, et al. 2010, ApJ, 715, L554 
Krichbaum, T. P, Witzel, A., Graham, D. A., et al. 1992, A&A, 260, L33 
Li, T.-P, & Ma, Y.-Q. 1983, ApJ, 272, L317 

Lombardi, S., Berger, K., Colin, P, et al. 2011, Proc. of 32nd ICRC, Beijin, 

China, arXiv:1109.6195 
Martin, P G., Thompson, I. B., Maza, J., & Angel, J. R. P 1983, ApJ, 266, L470 
Maza, J. 1979, PhD thesis, Univ. Toronto 

Moralejo, A., Gaug, M., Carmona, E., et al. 2009, Proc. of 31st ICRC, Lodz, 

Poland, arXiv:0907.0943 
Neronov, A., Semikoz, D. V., & Vovk, I. 2010, A&A, 519, L6 
Rolke, W. A., Lopez, A. M., & Conrad, J. 2005, Nuclear Instruments and 

Methods in Physics Research A, 551, 493 
Seyfert, C. K. 1943, ApJ, 97, L28S 

Vermeulen, R. C, Readhead, A. C. S., & Backer, D. C. 1994, ApJ, 430, L41 
Veron, R 1978, Nature, 272, 430 



IFAE, Edifici Cn., Campus UAB, E-08193 Bellaterra, Spain 

Universidad Complutense, E-28040 Madrid, Spain 

INAF National Institute for Astrophysics, 1-00136 Rome, Italy 

Universita di Siena, and INFN Pisa, 1-53100 Siena, Italy 

Technische Universitat Dortmund, D-4422I Dortmund, Germany 

Max-Planck-Institut fiir Physik, D-80805 Miinchen, Germany 

Universita di Padova and INFN, 1-35131 Padova, Italy 

Inst, de Astroffsica de Canarias, E-38200 La Laguna, Tenerife, 

Spain 

Depto. de Astroffsica, Universidad de La Laguna, E-38206 La 
Laguna, Spain 

University of Lodz, PL-90236 Lodz, Poland 

Deutsches Elektronen-Synchrotron (DESY), D- 15738 Zeuthen, 
Germany 

ETH Zurich, CH-8093 Zurich, Switzerland 
Universitat Wurzburg, D-97074 Wiirzburg, Germany 
Universitat de Barcelona (ICC/IEEC), E-08028 Barcelona, Spain 
Universita di Udine, and INFN Trieste, 1-33 100 Udine, Italy 
Institut de Ciencies de I'Espai (lEEC-CSIC), E-08193 Bellaterra, 
Spain 

Inst, de Astroffsica de Andalucfa (CSIC), E- 18080 Granada, Spain 
Croatian MAGIC Consortium, Rudjer Boskovic Institute, 
University of Rijeka and University of Split, HR- 10000 Zagreb, 
Croatia 

Universitat Autonoma de Barcelona, E-08193 Bellaterra, Spain 
Tuorla Observatory, University of Turku, FI-2I500 Piikkio, Finland 
Inst, for Nucl. Research and Nucl. Energy, BG-1784 Sofia, Bulgaria 
INAF/Osservatorio Astronomico and INFN, 1-34143 Trieste, Italy 
Universita dellTnsubria, Como, 1-22100 Como, Italy 
Universita di Pisa, and INFN Pisa, I-56I26 Pisa, Italy 
ICREA, E-080I0 Barcelona, Spain 

now at Ecole polytechnique federale de Lausanne (EPFL), 
Lausanne, Switzerland 
supported by INFN Padova 

now at: Centro de Investigaciones Energeticas, Medioambientales y 
Tecnologicas (CIEMAT), Madrid, Spain 

now at: Finnish Centre for Astronomy with ESO (FINCA), 
University of Turku, Finland 

HITS, Schloss-Wolfsbrunnenweg 33, 69118 Heidelberg, Germany 
UC Santa Barbara, CA 93106, Santa Barbara, USA 



4 



